PHYSICA  d 


ELSEVIER  Physica  C  368  (2002)  320-323  - 

www.elsevier.com/locate/physc 


A  multi- Josephson  junction  qubit 

S.P.  Yukon* 

Air  Force  Research  Laboratory,  Electromagnetics  Technology  Division,  Hanscom,  AFB  MA  01731,  USA 


Abstract 

We  have  designed  a  persistent  supercurrent  multi-Josephson  junction  (JJ)  qubit  whose  circuit  is  based  on  a  flattened 
JJ  triangular  prism.  The  Schrodinger  equation  for  the  ID  constrained  system  is  equivalent  to  the  Whittaker  Hill 
equation,  for  which  exact  solutions  have  been  found  [1].  Symmetric  or  antisymmetric  coupling  of  the  qubit  to  an  ex¬ 
ternal  magnetic  field,  will  excite  only  the  corresponding  symmetric  or  antisymmetric  terms  in  the  Hamiltonian.  This 
specificity  allows  coupling  to  a  system  bus  comprised  of  an  LC  resonant  loop.  We  indicate  how  separate  buses  might  be 
coupled  into  a  larger  branching  network.  Published  by  Elsevier  Science  B.V. 
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1.  Introduction 

The  qubit  circuit  is  based  on  a  discretized  and 
flattened  version  of  the  triangular  long  Josephson 
junction  (JJ)  prism  discussed  in  Refs.  [2-A]  and  is 
shown  in  Fig.  1.  The  long  JJ  prism  was  shown  to 
support  various  pairings  of  permanently  bound 
^o/3  kinks  (antikinks),  that  interpolate  between 
the  two  distinct  vacuum  states  of  the  system  at 
X  =  ±oo  (i.e.  in  the  two  JJ  end  triangles  of  the 
prism).  In  a  given  JJ  end  triangle,  the  two  vacuum 
states  correspond  to  plus  or  minus  circulating 
currents  J  =  Jc  sin(±7t/3)  around  the  triangle. 
Translated  to  the  flattened  prism  of  Fig.  1,  this 
implies  imposing  an  external  flux  of  =  ^o/3, 
where  the  flux  quantum  <Po  =  h/2e.  The  two  vac¬ 
uum  states  of  the  system  correspond  to  a  clockwise 
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circulating  current  in  the  left  (right)  triangle  and 
no  circulating  current  in  the  right  (left)  triangle.  In 
a  suitably  modified  system,  these  two  states  can  be 
taken  as  the  two  basis  states  for  a  qubit,  and  can 
be  shown  to  correspond  to  the  two  lowest  states 
((0|  ±  (l|)v/2  of  the  Hamiltonian  in  Eq.  (1). 

As  the  qubit  just  described  is  close  to  the  one 
that  has  been  investigated  by  Mooij  and  coworkers 
[5,6],  we  can  apply  the  results  of  their  investigation 
on  coupling  to  the  environment,  coupling  to 
SQUID  detectors,  and  the  low  inductance  limit  [7] 
to  the  system  proposed  here.  The  main  difference 
between  the  current  proposed  system  and  that  of 
Mooij  and  coworkers  is  the  spatial  mapping  of  the 
qubit  states  (which  in  turn  can  engender  different 
coupling  terms),  and  the  proposed  coupling  to  a 
Cirac  Zoller  type  bus  [8]. 

2.  Qubit  characteristics 

The  Hamiltonian  for  the  four  cell  qubit  of  Fig. 
1  is  given  by 
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Fig.  1.  Circuit  diagram  of  the  four  cell  flattened  JJ  triangular 
prism.  A  three  cell  variant,  with  one  central  cell  (corresponding 
to  the  opposite  prism  face  with  twice  the  external  flux),  is  also  a 
possible  qubit  circuit. 

H  =  2C{3{x<Po/2nf  +  {e^,/2nf) 

—  2itj{2cos(0)  cos((/)^/2)  cos((^^/2) 

X  (cos(x  +  0a/2)  -  2sin(0)  ?,m{(j)J2) 

X  sin(0^/2)  sinix  +  4>J2)  +  cos(2x) 

X  cos(0,  +  (1) 

where  we  have  defined  0  =  dj  —  0^/2  and  x  = 
lAi  +  (0,  -  0J/2  with  =  (0^1  ±  0^2)/2  and 
0^1 2  =  27t^^i_2/<Po,  along  with  the  diagonalizing 
plus  or  minus  combinations;  i/tj  =  (Ci  —  C2)/2, 
02  =  (^1  -  ^2)/2,  01  =  (Cl  +  C2)/2,  02  =  (Cl  +  C2)/ 
2,  as  well  as  0^^  =  (0^i  ±  0a2)/2  and  the  zero  in¬ 
ductance  limit  constraint  equations  {Ci  —  C2  —  C3  = 

0al>  Cl  ~  C2  ~  C3  =  0a2>  ~Cl  +  C2  =0^2>  C2  ~  Cl  = 

0^1},  where  0^,^,  =  27t<?,y,,,„/<?o  and  £'j  =  /c<?o/27i. 
The  junction  capacitance  C  and  critical  current  /c 
have  been  taken  to  be  equal  for  all  junctions.  We 
require  that  0^  =  0,  otherwise  mixing  of  even  and 
odd  X  rnatrix  elements  can  occur  upon  varying 
0„.  As  a  zeroth  order  Hartree  approximation 
we  assume  the  dynamical  variable  to  be  in  its 
ground  state,  and  make  the  replacement  cos(0)  ^ 
(0|  cos(0)|O)  =  ))  for  the  x  equation  of  motion, 
which  may  then  be  solved  exactly  [1].  An  harmonic 
expansion  of  the  0  potential  around  0  =  0  is  then 
made  with  cos(x)  ^  ((0|  cos(;t;)|0)  +  (1|  cos(x)|l))/ 
2  to  simplify  the  current  calculation.  Within  this 
approximation,  the  Hamiltonian  for  the  variable 
X,  may  be  written  as 

H  =  /2  -  2.Ej{2(0|  cos(0)|O)  cos(0,/2) 

X  (cos(x  +  <i>J2)  +  cos(2x)  cos(0J)},  (2) 


where  M  =  6(2C)(<Po/27t)^.  Replacing  p.^  =  Mx, 
the  momentum  conjugate  to  x,  with  —iJid/dx  then 
yields  the  Hamiltonian  operator  (in  units  of  2£’j), 

H/{2Ei)  =  -d^ldx^l2M^  -  {2i/ cos(0,/2) 

X  (cos(x  +  4>J2)  +  cos(2x)  cos(0,))}, 

(3) 

where  we  have  written  M^=  {2Ej)M  =  3{2Ej)/ 
(4(£c/2)),  with  Ec  =  eV2C. 

To  model  the  inclusion  of  rf  forcing  and  cou¬ 
plings  to  other  qubits,  we  expand  the  potential  in  a 
power  series  around  0°  =  0,  0°  0,  and  0°  ^  0, 

keeping  §0,-  terms  up  to  second  order.  The  0°  and 
§0,  magnetic  fields  can  be  supplied  by  small  loops 
or  by  control  lines:  0°  and  §0^  by  a  vertical  line 
bisecting  the  qubit,  0,.  by  a  horizontal  bisecting 
line,  and  0°  and  80,  by  the  dashed  line  in  Fig.  1, 
which  is  designed  to  maintain  0°j  2  =  0  and 
50^1 2  =  0.  In  the  Hilbert  subspace  spanned  by  the 
qubit  states  |0)  and  |1),  single  qubit  rotation  op¬ 
erators  that  are  first  order  in  the  50,  are  given  in 
terms  of  the  Pauli  spin  matrices  and  by: 

Z,  =  -))sin(0°)((O|cos(2x)|O) 

-  (l|cos(2x)|l))(50,,-F50^)ff„  (4a) 

=  -f/cos(0°/2)(O|sin(x)|l)50,,ff,,.  (4b) 

The  second-order  interaction  terms,  needed  for 
two  qubit  gates,  are  given  by 

A'b  =  -))sin(072)(O|  sin(x)|l)(50„50,/2)ff,„ 

(5a) 

A"b  =  -)7sin(0°/2)(l|  sin(x)|2)(50„50,/2)ff(, 

(5b) 

where  <t(.  =  (|2)(1| -F  |1)(2|).  Terms  proportional 
to  50J  will  contribute  to  a  small  shift  in  the  fre¬ 
quency  6{\Q^‘^<Dyc)  of  the  bus  oscillator  defined 
below. 


3.  The  control  Z  gate 

To  couple  the  qubits,  we  use  an  LC  resonant 
oscillator  as  a  system  bus,  based  on  the  trapped 


322 


S.P.  Yukon  /  Physica  C  368  (2002)  320-323 


ion  bus  scheme  of  Cirac  and  Zoller  [8,9]  (another 
scheme  based  on  an  LC  bus  using  a  direct  inter¬ 
action  has  been  described  in  Ref.  [10]).  The  bus  is 
coupled  to  each  qubit  with  a  gradiometer  type 
coupling  shown  in  Fig.  3,  that  nulls  or  §</)j 
single  qubit  excitations,  and  any  uniform  external 
perturbation  into  the  bus.  The  operators  for  flux 
and  charge  in  the  LC  oscillator  can  be  written  in 
terms  of  creation  and  annihilation  operators 
{a+,a}  as  [11]. 

d)  =  \/^Lffl/2(a  +  a+),  Q  =  i\/^/2Lffl(a  —  a+). 

(6) 

If  there  are  N  qubits  coupled  to  the  bus,  and 
a  fraction  fq<P  is  available  for  coupling  to  the 
qubits,  the  amount  of  flux  coupled  to  each  qubit 
will  yield  a  Jaynes  Cummings  type  [12]  in¬ 
teraction  terms  (5a)  and  (5b),  with  50,  =  2nfqX 
(^Lcu/2)'^^(a  -f  a+)/(A^^o)-  The  bus  frequency  ojlc 
is  chosen  so  that  neither  resonant  pulse  frequency 
(Y)a  =  fflio  -  fflLc  or  tt)b  =  ffl2i  -  CO-LC,  whcrc  noJij  is 
the  energy  needed  for  a  transition  from  level  j  to  i, 
is  close  to  a  qubit  or  bus  transition.  Solving  the 
Schrodinger  eigenvalue  equation  (3),  assuming 
€  =  0,  0°  =  0°,  r  =  -  cos(0.:)/cos(0°/2)  =  0.77, 
and  {2Ei)/h  =  200  GHz  yields  tt)io/27i  =  1.553 
GHz,  ffl2i/27i  =  15.09  GHz,  and  one  possible 
choice  of  coi^c/^ti  =11.5  GHz.  If  the  area  of  the 
bus  coupled  to  each  qubit  cell  is  taken  to  be 
A  =  I  [L^,  the  bus  inductance  for  «  1  can  be 
approximated  by  L  «  2Nl.25fiQ'/A  =  Nn  pH.  The 
amplitude  of  the  bus  oscillator  is  then 
~  1-67  X  IO^^^oVN.  For  r  =  0.77,  the  prod¬ 
uct  sin(0°)(l|  sin(;^)|2)/2  =  0.13.  Assuming  50,,  = 
27c(0.02<Po)/^o,  and  fq  =  1,  yields  a  total  X'^  am¬ 
plitude  of  1.7  X  /\/N  and  a  Rabi  frequency  of 
vr  =  34.3  MHz/AL2  for  the  slowest  (Jfg)  transition. 

4.  Sources  of  decoherence 

We  defer  a  detailed  discussion  of  decoherence, 
but  mention  that  in  addition  to  sources  discussed 
in  Refs.  [5,6],  the  bus  presents  a  new  source.  It 
can  act  as  a  linear  array  of  2N  small  1  loop 
antennas  with  an  RF  current  amplitude  of 
zrp  ss  <Pb/T  «  1.1  X  10“’  A  for  A  =  100,  yielding  a 
decoherence  time  from  radiation  losses  of  ^50  s. 


For  a  substrate  refractive  index  of  n  «  3,  the 
maximum  size  of  a  bus  path  is  ~Ai/10  ~  870  p, 
where  1,  =  l/n.  A  2^/20  square  branching  struc¬ 
ture  with  100  qubits  would  allow  the  qubits  to  be 
~48  p  apart.  While  the  large  RF  50„  excitation 
needed  for  two  qubit  gates  should  not  enter  the 
bus  due  to  geometric  cancellation,  its  mutual  in¬ 
ductance  effect  on  a  neighboring  (48  p  away)  qubit 
is  to  excite  a  gate  with  a  Rabi  frequency  of 
~500  Hz.  Reducing  N  to  50  and  orienting  alter¬ 
nate  qubits  to  be  mutually  perpendicular  could 
improve  this  by  a  factor  of  ^10“’. 


5.  Rabi  resonant  pulse  sequence  for  two  qnbit  gates 

The  complete  sequence  of  Rabi  resonant  pulses 
for  a  CZ  gate,  following  the  presentation  in  Ref. 
[9],  is  shown  in  Fig.  2.  The  notation  Xl^(ma) 
signifies  an  Ab  gate  pulse  of  n  frequency  coa  (de¬ 
fined  above)  acting  concomitantly  on  the  bus  and 
the  mth  qubit.  The  sequence  of  Ab  and  X'^  gates  in 
Fig.  2  has  been  successfully  modeled  using  the 
subset  of  24  basis  states  |m,«;7bus)  with  m  = 
{0, 1,2},  n  =  {0, 1,2,3},  y'bus  =  {0, 1}  and  the  sys¬ 
tem  parameters  discussed  above. 

As  the  number  of  qubits  on  a  given  bus  is  lim¬ 
ited,  we  consider  connecting  a  group  of  buses  in  an 
open  branching  network.  This  may  be  accom¬ 
plished,  as  shown  in  Fig.  3,  by  distinguishing  the 
qubit  that  connects  two  buses,  and  assigning  it 
the  sole  task  of  transferring  the  excitation  from  the 
first  bus  to  the  second  and  vice  versa.  To  carry  out 
a  CZ  gate  connecting  the  mth  qubit  on  bus  #1  with 
the  nth  qubit  on  bus  #2,  the  first  two  gates  of  Fig. 
2  are  carried  out  as  indicated.  The  next  gate 
"^B,/(®io  —  ®Lci)  transfers  the  excitation  on  bus  #1 
to  the  |1)  state  of  the  transfer  qubit.  A  second  Ag  , 
transfers  the  excitation  from  the  |1)  state  of  the 
transfer  qubit  to  bus  #2  with  excitation  frequency 
®io  —  fflLC2-  The  third  gate  in  Fig.  2  is  carried  out 
but  with  bus  #2  excitations  and  the  nth  qubit  on 
bus  #2.  The  final  step  transfers  the  bus  #2  exci¬ 
tation  back  to  the  original  bus  by  reversing  the  two 
original  X\^  gates.  The  sign  change  from  the 
round  trip  due  to  the  transfer  qubit  is  (— i)"^  =  1. 
Since  it  has  been  proved  [9]  that  the  combination 
of  one  qubit  gates  and  the  CZ  gate  is  sufficient  for 
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Fig.  2.  Gate  sequence  for  the  two  qubit  CZ  gate. 


Fig.  3.  Circuit  detail  for  the  coupling  of  two  resonant  LC  buses 
showing  a  transfer  qubit  (bi)  and  computational  qubits  (ai)  on 
bus  #1  and  (32)  on  bus  #2. 


producing  an  arbitrary  unitary  transform,  we  have 
demonstrated  that  the  coupled  bus  system  gates 
are  sufficient  for  a  branched  network  quantum 
computer. 
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